Effect of Inversion Symmetry on Incommensurate Order in Multiferroic RMn 2 05, 

R=rare earth. 



A. B. Harris[l], M. Kenzelmann[2], Amnon Aharony[3], and 0. Entin-Wohlman[3] 
[1] Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104 
[2]Laboratory for Solid State Physics, ETH Zurich, 
CH-8093 Zurich, Switzerland and Laboratory for Neutron Scattering, 
ETHZ & PSI, CH-5232 Villigen PSI, Switzerland and 
[3] Department of Physics and the Use Katz Center for Meso-and Nano-Scale Science and Technology, 
Ben Gurion University, Beer Sheva 84105 ISRAEL 
(Dated: March 7, 2008) 

Starting from the irreducible representations of the group of the wave vector we construct the spin 
wave functions consistent with inversion symmetry, neglected in the usual representation analysis. 
We obtain the relation between the basis functions of different members of the star of the wave 
vector. We introduce order parameters and determine their transformation properties under the 
operations of the space group of the paramagnetic crystal. The results are applied to construct 
terms in the magnetoelectric interaction which are quadratic and quartic in the magnetic order 
parameters. The higher order magnetoelectric interactions can in principle induce components of 
the spontaneous polarization which are not allowed by the lowest order magnetoelectric interaction. 
We also obtain the relation between the spin wave functions of the incommensurate phase and those 
of the commensurate phase which lead to analogous relations between the order parameters of these 
two phases. 

PACS numbers: 75.25.+z,75.10. Jm,75.40.Gb 



I. INTRODUCTION 



The problem of determining the symmetry of incom- 
mensurate (IC) magnetic order from diffraction exper- 
iments is an old one and is the subject of several well- 
known reviews i 2 - The reviews are based on the idea that 
the spin structure that develops at a continuous transi- 
tion must transform like an irreducible representation (ir- 
rep) of the group of operations which leaves the IC wave 
vector q invariant. 3 However, perhaps surprisingly, these 
standard references do not exploit additional restrictions 
that are due to inversion symmetry when that opera- 
tion is not a member of the group of the wave vector. 
Although the group theoretical formalism for doing this 
has been described^ and these restrictions had previ- 
ously been used to aid in structure determinations^^^ 
the effect of inversion symmetry is often not included in 
the classification of possible magnetic structures. 

Here we perform the requisite analysis for the star of 
wave vectors of the IC phase a 10 ' 11 ! 12 of the "125" sys- 
tems, RMn2 05, where R is a rare earth ion, which may 
be magnetic or not (e. g. when R is yttrium). The inter- 
est in these materials stems from the fact that they ex- 
hibit ferroelectricit y 13 i 14 ' 15 i 16 whose onset coincides with 
a magnetic ordering transition! 17 i 18 i 19 i 20 i 21 We show that 
when inversion symmetry is taken into account, there are 
about half as many degrees of freedom that describe the 
basis functions of the irreducible representations com- 
pared to an analysis when inversion symmetry is over- 
looked. Even when an unrestricted fit (not taking ac- 
count of any symmetry) is performed, 22 it is useful to 
have the results of the present paper to see if the hy- 
pothesis of a single irrep 2 ^ holds. Thus, it is clear that 



magnetic structure determination using an approach that 
includes inversion symmetry will lead to an increase in 
the accuracy of the structure determinations. Finally, 
this approach leads naturally to the introduction of or- 
der parameters which have symmetry properties that we 
explicitly display and in terms of which a Landau ex- 
pansion was developed for a number of systematise and 
which has led to a generic magnetoelectric (ME) phase 
diagram for the 125'sj 2 ^ The purpose of the present paper 
is to a) analyze the symmetry of the various IC phases, b) 
show how the symmetry implies relations between order 
parameters of different symmetry magnetic phases, and 
c) analyze the symmetry of the ME interactions which 
explain the appearance of ferroelectric order at some of 
the magnetic phase transitions. 

Briefly, this paper is organized as follows. In Sec. II 
we list the results obtained using the canned program 
MODY for the IC phase and we show how to modify 
this to take account of inversion symmetry. Here order 
parameters are introduced as the complex amplitudes of 
the spin wave functions. In Sec. Ill we show how, hav- 
ing obtained the basis functions for one member of the 
star of q, one can determine the basis functions for all 
the other wave vectors in the star of q. Here we also 
determine how the order parameters transform under all 
the operations of the space group. Having determined the 
symmetry properties of the order parameters we are able, 
in Sec. IV, to construct the lowest order (trilinear) ME 
interaction which explains the orientation of the observed 
magnetically induced spontaneous polarization. Here we 
show that higher order and Umklapp ME interactions can 
lead to small contributions to all components of the spon- 
taneous polarization. In Sec. V we discuss how the basis 
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FIG. 1: (a): The ME phase diagram of ErMn 2 5 .— Here 
(X,0, Z) indicates the nature of the wave vector. If X = C 
(Z = C), then q x = 1/2 (q z = 1/4). If X = I (Z = I), then 
q x (qz) is IC, but close to 1/2 (1/4). The dashed lines indi- 
cate temperatures at which an anomaly in the &-component 
of the dielectric constant was observed. P\\b indicates that 
the system has a spontaneous polarization aligned along b 
(for T < 39 K). (b): Same for HoMn 2 5 ^ For T > 39 K, 
q z < 1/4 and for T < 20 K, the (1,0,7) phase has q z > 1/4 
and the system is either paraelectric or weakly ferroelectric. 



functions in the IC phase with q x ^ 1/2 connect to those 
in the adjacent q x = 1/2 phase. Here we also analyze 
the symmetry of the special multicritical point for which 
q x = 1/2. In Sec. VI we briefly summarize the results of 
this paper. 

II. CALCULATION 
A. Results without Inversion Symmetry 

The lattice structure of the 125's was determined by 
Quezel-Abrunaz et al. 25 to be that of the orthorhombic 
space group Pbam (#55 in Ref. [26h . In Table|T]we list the 
general positions in the primitive unit cell which define 
the symmetry operations of the space group Pbam and 
in Table I I 27 ' 28 i 29 we give the actual positions of the ions 
for the 125 systems. 

The magnetic and dielectric phases occuring in the 
125's are more complicated and we give a brief overview 
of them here. In Fig [Iji and [lb we show the ME 
phase diagrams of ErMr^Os (taken from Ref. [l7l ) and 
HoMn205 (taken from Ref. [la ) which exhibit the si- 
multaneous ferroelectric and magnetic phase transitions. 
When cooled from the paramagnetic phase, the 125's de- 



TABLE I: General positions within the unit cell for space 
group Pbam expressed as fractions of the orthorhombic lattice 
constants.— This table defines the space group operations on 
r = (x,y,z). Here 2 a is a two-fold rotation (or screw) about 
the a axis and m a p is a mirror (or glide) a/3 plane. 
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TABLE II: Position t„ (in units of lattice constants) of the 
nth magnetic ion in the unit cell. (These values are for 
HoMn2 05 , 27,28 but are approximately the same for the other 
125's.^) Sites 1-4 are for Mn 3+ , 5-8 are for Mn 4+ and 9-12 
are for R 3+ ions. 
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velop IC order at about 45 K in a paraelectric phase 
described by the wave vectors whose star consists of 
q± = [(1/2 — 8, 0, ±(1/4 + e)] and their negatives, where 
6 and e are of order 0.05 or l ess ii^^^^^^^^ 
in reciprocal lattice units (rlu's). Upon further cooling 
some 125's, such as ErM^Os (shown in Fig. la) ) 17 i 21 
YMnzOf^^&SO and TmMnzOs,-^^ exhibit a ferro- 
electric (1,0, C) phase in which e = 0, before entering a 
(C, 0, C) phase in which 8 — e — 0. Other 125's, such as 
TbMn 2 5f ii^^HoMn 2 05 (shown in Fig. lb)^^^^ 
and DyMn 2 5 ^i2i^i26 g0 directly from the (7,0,7) 
phase into the (C, 0,C) phase without the appearance 
of the (7, 0, C) phase. At lower temperature the 125's 
follow various scenarios in which the magnetic structures 
may be either IC or CM with a long period and they may 
or may not be ferroelectric. For a review of the properties 
and Landau theory for 125's see Ref. 1371 

Here we give a symmetry analysis of the allowed mag- 
netic structures in the (7,0,7) or (7,0, C) phases. (A 
detailed symmetry analysis applicable to the phase with 
q x = 1/3^ indicated that this phase was described by a 
two dimensional (2D) irrep and therefore could be char- 
acterized by two complex- valued order parameters^ we 
will call o i and (72 ■) The symmetry of the phase when 
q x 7^ 1/2 is different. The group of this wave vector 
contains unity E and the glide m ac which leaves the b- 
component of the wave vector invariant. Thus we have 
two one dimensional (ID) irreps, which we label T e and 
r o ("e" for even and "o" for odd). In particular, since 
the star of the wave vector contains four vectors, ordering 
within each irrep is described by four complex-valued or- 
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TABLE III: Symmetry-adapted basis functions for wave 
vectop 3 ^ q+ = (q x ,0,q z ) which transform according to the 
irreps T e and r o , where A = exp(— niq x ). We have not yet 
included the effect of inversion symmetry. 
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der parameters.— The allowable wave functions are the 
basis functions of the irreps which transform appropri- 
ately. These basis functions are actually eigenvectors of 
m ac with eigenvalues +A* (for T e ) and —A* (for r o ), 
where A = exp(— iitq x ). Since each irrep is contained 
18 times in the original reducible representation gener- 
ated by the three spin components of the 12 magnetic 
sites in the unit cell (here we assume that R is mag- 
netic), each wave function contains 18 independent free 
complex-valued parameters. These wave functions are 
listed in Table HIB ^ and they are in agreement with (i. 
e. are a reparameterization of) the results of the MODY 
program. 

To illustrate the transformation laws, we perform a 
partial check that the vectors in Table IIIII are indeed 
eigenfunctions of m ac . Note that we use so-called "unit- 
cell" Fourier transforms whereby 40 



S(R,n) = S(q,n)e- 2,r4q - R + c. c. 



(1) 



where n labels the sublattice and R locates the unit cell. 
A transformation O takes the "initial" basis function into 
a "final" basis function. If a prime indicates "final," i. e. 
"after transformation," then S"(R/, n/) denotes the spin 
of sublattice nj in the unit cell at R/ after transforma- 
tion. This quantity is obtained by applying the transfor- 
mation to the spin at the initial location R^ + T ni . Thus 
for transformation by m ac we write 



S'(q,n f ) = £ Q S Q (q,n;)e 



27riq-[R / -R i ; 



(2) 



where £ Q is the appropriate factor for the mirror oper- 
ation m ac on the components of a pseudovector: t; y — 
—£ x = —£, z = 1. We will check that the basis vector of 
irrep r e is an eigenfunction of m ac . Note that under m ac 
when the initial sublattice index is = In — 1, then the 



final sublattice index is rif 



In and vice versa. Thus 
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Thus ip(T e ) is an eigenvector of m ac with eigenvalue A*. 
In the other irrep, the fact that A is everywhere replaced 
by —A ensures that tp(T ) is an eigenvector of m ac with 
eigenvalue —A*. 



B. Effect of Inversion Symmetry 

Now we modify the above results to take account of 
inversion symmetry. A straightforward, if clumsy, way 
to do this is to use the fact that the inverse susceptibil- 
ity matrix becomes singular at a continuous phase tran- 
sition, which implies that one of its eigenvalues passes 
through zero. We wish to see what restrictions symme- 
try places on the associated critical eigenvector. We write 
the quadratic terms in the free energy F 2 in the form 



Fo 



1 



(4) 



where T is the inverse suscptibility matrix. Instead of 
considering the quadratic form in the original spin vari- 
ables, we consider the quadratic form in terms of the vari- 
ables of Table Hill So the matrix T is an 18 dimensional 
Hermitian matrix operating on an 18-component vector 



4 



W(r) which we write as (si,ti,S2,t2,s 3 ,t 3 ), where the 
s's and t's are three component subvectors taken from 
Table InH Thus 



(5) 



Because the paramagnetic phase has symmetry under 
spatial inversion Z, we must have^^. 



(6) 



f 2 = -pra]tjF[:ra] = 



for all values of the spin coordinates. 

To implement this we note that for transformation un- 
der X, the result follows a logic similar to that leading to 
Eq. (J21), namely^ 



(7) 



where again the prime indicates the value after transfor- 
mation by X. Note that inversion relates sites (1,4), (2,3), 
(5,8), (6,7), (9,12), and (10,11). Now use Eq. © to get 
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These simple results arise because we reparametrized 
with an eye to avoiding complexity. 

The eigenvalue equation for the 18x18 matrix T can 
be represented as 



(10) 



where each entry of the matrix is itself a 3 x 3 submatrix. 
Now we identify the symmetry of this matrix imposed by 
inversion, via Eq. We have 
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G jit j\til , 

which implies that Ay = Gji, so that G = A = 
since A is Hermitian. Similarly, P — L* and U 
Consider 

Bij s iitji — Bij[Xsn] [Xtji] = B i jtnSj 1 



(11) 

A*, 
= S*. 



(12) 



which implies that Bi 



Bji. Thus Bt 



B*. Likewise 



Mt = M* and Tt = T*. Furthermore 

GijS il Sj2 = Cij[Xsn] [XSj2\ = Cijtntj2 

— ft ../* 

— j2 l ± 

which implies that I*j = C'ij. Also 

E ij s *i s j3 = Eij[Xsa]*[Lsj3,] 



(13) 



(14) 



which implies that K* 
Ne 2 ™K Also 



Ee 2 ™«*. Similarly R* 



DijS^tj-2 = D lj [Xsii]*[Xtj 2 ] = DijtiiSj 2 



(15) 



so H* = D. Also 



Jijt^Sjs — Jij[Xta]*[Xsj 3 ] 
= JijSnt* 3 [e 
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(16) 



Similarly O* 



which implies that F* 
Qe 2 ™ 9 * . 

Using all these relations we see that the matrix J- must 
be of the form 



ABC 

B* A* D* 
C* D L 
D* C M* 
E* J* N* 



D E J*A* 
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M N Q*A* 

L* Q N*A* 

Q* S T 



JA EA QA NA T* 



(17) 



where A = exp(27rig z ). Now consider this matrix oper- 
ating on a vector of the form 



* = [ P ,p*,ii>,ip*,xAx*] 



(18) 



One can show that T^S is a vector of the same form as 
. This means that any eigenvector can be taken to be 
of this form and the eigenvalue equations are 

EV + JV + Nty + QV + S X + TAX* = Ax 
C f p + Dp* + Lif> + MtP* +N X + Q*X* = AV 
Ap + Bp* + CV + DV* + Ex + J*X* = Xp ■ (19) 

(The other three equations are the complex conjugates 
of these.) These give rise to 18 simultaneous equations 
for the real and imaginary parts of the three-component 
vectors p, and X - 

The point is that the permissible form for an 18- 
component eigenvector is restricted by inversion sym- 
metry. The critical eigenvector is the one whose eigen- 
value first passes through zero as the temperature is low- 
ered. As the temperature is further lowered, we may have 
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a small amount of admixing of noncritical eigenvectors 
into the critical eigenvector due to higher than quadratic 
terms in the free energy. However, these admixtures will 
only be within the same irrep unless one crosses a phase 
boundary. 

Since the eigenvalue problem is in a complex vector 
space we write critical eigenvector as 

* = e ^[ P) p*,T/,,^*,x,A X *] , (20) 

where the phase <f> is arbitrary (as far as the quadratic 
terms are concerned) and the other Greek letters are 
three component vectors. In Tables IIVI and [V] we tab- 
ulate the results. In so doing we have introduced the 
complex- valued order parameters cr(T), such that 

a(T) = |cr(r)|e^( r) . (21) 

To avoid overparametrizing we specify the normalization 

4 J2 E Kn\ 2 = l- (22) 

a— x,y,z n— 1,2,3 

Including inversion symmetry we have 9 complex-valued 
s parameters and one complex valued order parameter 
cr e (q+), so that we have 19 real valued parameters (tak- 
ing account of the normalization of the s's.), whereas 
without taking account of inversion symmetry we would 
have had 36 real valued parameters to determine from a 
fit to diffraction data. 

One may notice that we could have said that the 18- 
component eigenvector of u's was of the form 

$ = e^[n,-7T*,T,-T*^,-H*} (23) 

and indeed the eigenvector is equivalent to this form be- 
cause if you multiply the previous eigenvector 'J/ by i, it 
will be exactly of the form of <&. 

The comparison with Ni 3 V 2 8 6 < 8 < 9 (NVO) and 
TbMnO^ (TMO) is significant. In the case of NVO 
the magnetic Ni sites are of two types, spine and cross- 
tie£ All sites of the same type are related to one another 
by a symmetry operation which leaves the wave vector 
invariant. It happens that the Wyckoff orbit of this set 
of operators generates the entire set of spine sites and 
also separately the entire set of cross-tie sites. In that 
case inversion (which does not leave the wave vector in- 
variant) fixes all the relative phases^^ (The phases are 
not necessarily the same, but they are fixed.) In the case 
of TMO the Mn sites form a Wyckoff orbit of the sym- 
metry operations that leave the wave vector invariant, 
but the Tb sites break into two orbits. In this case in- 
version fixes the relative phases within the Mn orbit and 
within a single Tb orbit. Inversion connects the two Tb 
orbits. As a result the amplitudes of the two Tb orbits 
are fixed to be the same and they have phases which are 
the negatives of one another, but the magnitude of this 
phase is arbitrary^ 9 - Here the Mn 3+ , Mn 4+ , and RE sites 
each break up into two orbits which are interconnected 
by inversion. So it is not surprising that this situation is 



TABLE IV: Symmetry-adapted spin wave functions for wave 
vector q+ = (q x ,0, q z ) which transform according to the irrep 
r e , where A = exp(— niq x ), A — exp(2-7rig z ), and <r e is the 
complex-valued order parameter. We require the normaliza- 
tion of Eq. (|22[) . Otherwise, all constants assume arbitrary 
complex values. Here we include the effect of inversion sym- 
metry. 







S(q,l) = 


&e(q.+)(s x i, Syi, s 2i ) 


S(q,2) = 


— <T e (q+)\(Sxl, —Syl, Szl) 


S(q,3) = 


-<T e (q+)A*(s* lj -s£ 1 ,si 1 ) 


S(q,4) = 


o"e(q+)(sJi,s* 1 ,s* 1 ) 


S(q,5) = 


0- e (q + )(s l2 , S y 2, Sz2) 


S(q,6) = 


-cr e (q + )A*(s a; 2, — s y2 , s z2 ) 


S(q,7) = 


-o- e (q + )A*(s* 2 , -s* 2l s* 2 ) 


S(q,8) = 


o- e (q+)(s* 2 , s*2, s* 2 ) 


S(q,9) = 


o"e(q+)(srf, s v3, s z s) 


S(q, 10) = 


-cr e (q + )A(s a; 3, -s y3 , s z3 ) 


S(q,ll) = 


-<j e (q + )A*A« 3 , -s* 3 , s* 3 ) 


S(q, 12) = 


<7 e (q + )A(s*3,S*3,S*3) 



TABLE V: As Table HV] but for the irrep E and we require 
the normalization of Eq. ([22} with s replaced by u. 





^(Lo) 


S(q, 1) 




ffo(q+)(« a! i, 


Uyl, U z l) 


S(q,2) 




a (ci + )\(u x i, 


—U X 2, Uzl) 


S(q,3) 




cr (q + )A*« 1 


-u yl ,u zl ) 


S(q,4) 




cr (q+)« 1 , 


u "yl, u zl) 


S(q,5) 




<r (q.+)(u m 2, 


Uy2,U z 2) 


S(q,6) 




cr (q+)A* (u x 2 


— U y 2, Uz2) 


S(q, 7) 




cr (q + )A*«2 


-U y 2,U*z2) 


S(q,8) 




c (q+)(u* 2 , 




S(q,9) 




cr (q + )(u a , 3 , 


Uy3,Uz3) 


S(q, 10) 




a (q+)\(u x 3, 


— Uy3,Uz3) 


S(q,ll) 




cr (q + )A*A(M* 3 


, -U* 3 ,wJ 3 ) 


S(q,12) 




a a (q+)A« 3 


! u y3i u *z3) 



like that of the Tb sites in TMO: the magnitudes of the 
two related orbits, which according to MODY were un- 
related, are now, by virtue of inversion symmetry, fixed 
to be the same. 



III. DISCUSSION 

A. Order Parameters 

It is natural to introduce order parameters because as 
the temperature is reduced into the ordered phase, the 
critical eigenvector is nearly temperature-independent 
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so that 



Similarly 



except for a change in its normalization, governed by Similarly 
the magnitude of the order parameter. Furthermore, the 
phase of the complex order parameter is either a free 
variable or, if it is fixed, it is only fixed by subtle effects 
of higher-than-quadratic terms in the free energy. So the 
order parameter describes properly the low energy sector 
of the free energy. 

Note that our definition of the order parameter is such 
that if one is given the spin wave function over all the 
sublattices it is possible to uniquely determine both the 
phase and the magnitude of the order parameter, except 
that it could be multiplied by —1. (But that indetermi- 
nacy is inherent for this order parameter symmetry.) To 
make this unique identification from a knowledge of the 
wave functions, the wave functions must be first put into 
the canonical form of Tables HVl and fVl In so doing, the 
normalization condition has to be obeyed. Then the pref- 
actor will be the desired order parameter. Note that the 
phase is fixed by having the first and fourth components 
written in terms of complex conjugates. This type of 
identification would not be possible for a one-component 
complex variable. 

It should be noted that the order parameter inherits 
the symmetry of the full wave function. Having the basis 
functions for— q + = (q x ,0, q z ) we now obtain the basis 
functions for the other wave vectors in the star of q. We 
first obtain the basis functions for q = (— q x , 0, q z ) for 
irrep r e . The most general basis function for irrep F e for 
this wave vector will be of the form of Table IIVI with q x 
replaced by — q x , i. e. with A replaced by A* and, for 
notational convenience, s Q)Tl replaced by t a , n , However, 
this is not the basis function we want. We want the 
particular basis function which is obtained from that of 
q+ by a symmetry operation which takes q+ = (q x , 0, q z ) 
into — q = (—q x , 0, q z )r^ because it is this basis function 
which results from the actual interaction between spins. 
In other words, we want to relate t a n to s a n . To do 
this we now study the transformation of the spin Fourier 
transforms. so t h a t 

We first consider transformation by 2 C which takes 
q = (q x ,0,q z ) into q' = (-q x , 0, q z ) = -q_, where here 
and below we use a prime to indicate a quantity after 
transformation. We have that Similarly 



S' a (R f A) = p a S a (R l ,l) 



(27) 



with r, = (X,Y,Z) + n, 17 = (X - 1, Y - 1, Z) + r 4 
which gives 



U* -2m[(-q x ,0,q z )-(-X-l ,-Y-l,Z)] 
n l al e 



p a crs a ie 



-2m[(q x ,0,q^}-{X,Y,Z)] 



I \ 2 1* 
r] A t al = p a (TS a i 



S' a (Rf, 5) = p a S a (R 2l 5) 



(28) 



(29) 



(30) 



with r, = (X, Y, Z) + 1-5, r f = (X — 1, Y, Z) + r 5 which 
is 



v > ta2e -^ii(-q*,o,q z )-(-x-i,-Y-i,z)} 



p a crs a2 e 



-2iri[(q x fl,q,)-(X,Y,Z)] 



so that 



Similarly 



t/A 2 ^ = p a os a2 



Saffi-fl 7) — PaSailli, 7) 



(31) 
(32) 
(33) 



with r< = (X, Y, Z) + tt, T f = (X, Y - 1, Z) + t 7 . In 
using Table IIVI we must replace A by A* to convert the 
table for the wave vector q'. Thus Eq. (f3"3")) yields 

= p a Z a \*as* a2 e- 2 ™^*> ^< x ' Y < z V (34) 



(35) 



S' a (R f ,i) =p a S a (RiA) 



(24) 



S , ;(R / ,9)=p Q 5 a (R i ,12) 



(36) 



where p x = p y = —p z = — 1. We now write this in 
terms of Fourier components using Eq. (|TJ) . The initial 
position is = (X, Y, Z) + T4 and the final position 
is Vf = (X — 1,Y — 1,Z) + n which gives [with rj = 
CT e(-q) = cr e (q_)* and a = cr e (q+)] 



with r; = (X,Y,Z) + T 12 , r f = (X-1,Y-1,Z) + T9, 

which gives 

^ a3e -2irt[(-g*,0,g«)-(--X--l,-y-X,2)] 
= p a aAs* a3 e- 2 ™^*> >^< x > Y > z V (37) 



Tj r taie -2ni{(-<l*,0,q*H-X-l,-Y-l,Z)] 
n as* p-2irt[(?»,0, to )-(X ) Y ) Z)] 



So with exp[— 2wiq x ] = A 2 , we have 
■q'\ 2 t a i = p a crs* al . 



(25) 



(26) 



so that 



Similarly 



3 • 



S' a {K f ,l2) = p a S a (IL i ,9) 



(38) 
(39) 
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TABLE VI: Amplitudes of the basis functions for the irrep 
r„ for the star of q, where p a = ( — 1, — 1, 1). Here we give 
the basis functions for sublattices 1, 5, and 9. The remaining 
amplitudes are found by the appropriate modification of Table 
IIVI for the wave vector in question. For the irrep F , replace 
all the s's by it's and the remaining amplitudes are found by 
the appropriate modification of Table IVl for the wave vector 
in question. 



q + = (qx,0,q s ) 


-Sal S a 2 Sa3 


q = {-qx,o,q z ) 


PaS* al p a S a 2 Ap a S a3 


q = (qx,0, ~qz) 


PaSal paS a2 A p a Sa3 


-q + = (-<7z,0, -q t ) 


s al s ct2 s a3 



the sign of the wave vector leads to complex conjugation 
of the eigenvalue. 

We consider next the effect of 2 C on the order param- 
eters. In Eq. (|43|) we found that under 2 C the new value 
of vq is X* 2 v\. Since the prefactor A* 2 does not depend 
on q z and it was obtained without specifying the irrep, 
we see that the prefactors in the last four columns of the 
second row are the same. The prefactors of the first four 
entries of this row are obtained from the last four entries 
by complex conjugation. 

Next we consider the effect of inversion on the order pa- 
rameters. This discussion is simplified by having in hand 
the results of Table IVII Note that X does not change the 
orientation of the spin, because spin is a pseudovector. 
So under X we have 



with = (X,Y,Z) + t 9 , T f = (X- l,y-l,Z) + Ti2, 
which gives 

A^'£* 3e -27ri[(-fe,0,g*H-X-l,-y-l,.Z)] 

-2iri[(q x ,0,q z )-{X,Y,Z)] uq\ 



— Pa&Saze 



so that 



r)'\ 2 At* a3 = p a Aas a3 



(41) 



Equations {26 ]) . rf29 j) . (|32 j f35jl . (j38 L and (jlB yield 



tal — PaS* a i , t a 2 — p a S a 2 , t a3 — Ap Q S* 3 , (42) 



and 



V 



X* 2 a 



(43) 



There is an equivalent solution in which all the trans- 
formed quantities are multiplied by —1. This ambiguity 
is unavoidable because it is inherent in the symmetry of 
the order parameter. Using Eq. (|42[) and the fact that 
the basis functions for — q are the complex conjugates of 
those for q we obtain the results of Table IVII The rela- 
tions for the basis functions of irrep T a are the same as 
for r e , so Table [VTI also applies for r o . 

We now obtain the transformation properties of the 
order parameter under all the symmetry operations of 
the space group (except translations). For this discussion 
it is convenient to introduce an order parameter vector 
v whose components are the various order parameters: 

vi = cr e (q + ) , v 2 = CT e (q_) , v 3 = cr (q + ) , 
v 4 = CTo(q-) , v 5 = cr e (-q+) , Vq = cr e (-q_) , 
V7 = cr c (-q+) , i> 8 = cr (-q_) . (44) 

The transformation properties of the vector v are given 
in Table IVII1 whose construction we now discuss. The 
row of m ac is obtained by using the fact that the basis 
vector of irrep T e for wave vector q = (q x ,0 7 q z ) is an 
eigenvector of m ac with eigenvalue A*. The eigenvalue 
for irreps T e and r o have opposite signs, and changing 



(45) 



where r 4 = (X, Y, Z) + r 4 and t } = (X - 1, Y - 1, Z 
1) + Ti, which gives [with r) = er e (— q) and a — <J e (q)] 

/ / -2m[{-q x fi,-q z )-{-X-l,-Y-X,-Z-l}} 
V s al e 

i '%'* p 2m[(-q x ,0 ,- q ,).(-X-l,-Y-l ,-Z-l)] 
Til * Q iC 



crs^e 



+cr s a \e 



-2m[(q m ,0,q z )-(X,Y,Z)] 
2ni[{q x fl,q z y{X,Y,Z)] 



(46) 



This has to be an equality for all integer X, Y, and Z. 
Also s' al = s*j (from Table [*VT]) . so we find that 



r/A 2 A* 



(47) 



Thus A* Avi is the entry under V5 in the third row. Hav- 
ing this result, one can construct the other entries in this 
row by noting the dependence on q x and q z . 

The other rows of Table IVIII arc found by using the 
multiplicative properties 



m ab = 2 C X , 
m bc = 2 a X , 



2a ^ac^ab ? 

2 fc = m ac X . (48) 



IV. MAGNETOELECTRIC INTERACTION 

Now we discuss the form of the ME coupling in the 
phases with q x ^ 1/2, i. e. in the (I, 0,1) and (1,0, C) 
phases. In the first subsection we will discuss the tri- 
linear ME interaction which involves the lowest num- 
ber (two) of magnetic order parameters. In succeeding 
subsections we will discuss higher order ME interactions 
which involve a product of four magnetic order parame- 
ters. These higher order terms yield components of the 
spontaneous polarization which are allowed by symmetry 
but are not present in the trilinear interaction. However, 
these higher-order terms are probably small for two rea- 
sons. Firstly, in the IC phases which occur at high tem- 
peratures near the paramagnetic phase, the order param- 
eters are small. Secondly, most microscopic models of the 
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TABLE VII: The first column gives the operation O and the column headed v n gives the result of Ov„ where v is given in Eq. 
(|44|) . The last column gives the eigenvalue of dV lvA /dP b in Eq. (|53p under the operation O. 



o 


Vl 


V2 


«3 


V4 


«5 


^6 


«7 


^8 


dV int /dP b 


m ac 


X*Vi 


X*V 2 


-A*v 3 


— X*V4 


Xv s 


Xv 6 


— A«7 


-Au 8 


-1 


2 C 


X 2 v 6 


X 2 v 5 


A 2 « 8 


X 2 v 7 


X* 2 V 2 


A*V 


A*V 


•> *2 
A V3 


-1 


T 


A 2 A*« 5 


X 2 Av 6 


A 2 A*-u 7 


X 2 Av 8 


A* 2 A«i 


X* 2 A*V2 


A* 2 A« 3 


A* 2 A*t;4 


-1 


m bc 


\v 6 


A« 5 


-Xv$ 


— Xvt 


X*v 2 


A*vi 


— A*t>4 


-A*n 3 


+1 


2a 


\*A*V 2 


A*A«i 


-X*A*v 4 


-X*Av 3 


AAv 6 


AA*v 5 


-AAi; 8 


-AA*u 7 


-1 


m ab 


A*v 2 


Avi 


A*v 4 


Av 3 


Av 6 


A*« 5 


A« 8 


A*w 7 


+1 


2t 


XA*v 5 


XAv a 


-AA*v 7 


-A Au g 


X*Av 1 


A*A*« 2 


-A*Au 3 


— A*A*« 4 


+1 



ME interactio n 30 i 42 i 43 i 44 i 45 treat (within lowest order per- 
turbation theory) a trilinear Hamiltonian involving two 
spin variables and one displacement variable. However, 
to obtain these higher order phenomenological interac- 
tions probably involve processes of higher order in some 
small parameter such as t/U or X/U, where t is a hop- 
ping matrix element, A is the spin-orbit constant, and U 
is a Coulomb interaction. 



A. Trilinear ME coupling 

Initially we will consider the lowest order (trilinear) 
ME coupling. We start by considering the case when only 
the wave vectors ±q + = ±(1/2 — 5, 0, 1/4 + e), where e 
may or may not be zero, are involved. The interaction 
of lowest order in the magnetic order parameters which 
conserves wave vector and is time-reversal invariant is of 
the form&£i£ 

Vint = ^2 C 7afc cr a (q + )cr 6 (-q + )P 7 , (49) 
7a, b 

where a and b assume the values "e" and "o", P is 
the spontaneous electric polarization and 7 labels the 
component^ Using Table IVIIi one sees that terms in 
Vint with a = b are not allowed by inversion invariance. 
If one has only a single irrep present, then one can always 
redefine the location of the origin so as to have inversion 
symmetry with respect to that new origin and hence such 
a phase can not exhibit magnetically induced ferroelec- 
tricity. If both irreps are present, then we write 

Vint = ^[c 7 cr e (q + )cr (-q+) + c 7 cr e (-q + )cr (q + )]P 7 , 

7 

(50) 

and inversion invariance forces c 7 to be pure imaginary: 
c 7 = ir~ n where r 7 is real. Then 

Vint = i^r 7 [cr e (q + )cr (q + )* - fT e (q+)V (q + )]P 7 . 

7 

(51) 



From Table IVlTl one sees that the square bracket in this 
equation changes sign under m QC , so P 7 must also change 
sign under m ac in order for V; n t to be invariant under m ac . 
Thus c 7 can be nonzero only for 7 = b, as is observed. 
If we set or(q+) = |or(q+)| exp(i(/)r), then we have the 
result 

Vint = 2rsin(0 o - e )Pf,|cr e (q + )f7 o (q + )| . (52) 

However, this is not the whole story because we must 
include the terms involving the other wave vectors in 
the star of q. (Indeed it is possible that in the highest 
temperature paraeletric IC phase there is a simultaneous 
condensation of the order order parameters of both wave 
vectors q± 1^) Since we have already incorporated the ef- 
fect of X and m oc , it only remains to use 2 C to obtain the 
other terms which make up the invariant interaction. To 
do that we use the results given in Table IVIII which give 
2 c Cn(q+) = A 2 er„(q_)*, for n = o or e, and, of course, 
2 c p, = — Pfe. Thereby we obtain the complete result for 
Vint: 

VI n t = ir^2{(T e (q v )<j (q. n )* - <j e (q v )*<j (q v )}P b . 

77=± 

(53) 

At this order one needs the simultaneous presence of both 
the e and o irreps to have ferroelectricity. (However, be- 
low we find that a polarization along c can be induced 
by Umklapp ME interactions by a single irrep. But this 
scenario is unlikely. 24 ) Note that from this interaction 
the spontaneous polarization P is aligned along the b 
axis irrespective of which wave vector condenses. How- 
ever, the sign of P depends on how the signs of the order 
parameters are chosen (i. e. how symmetry is broken) 
when a and/or a e order. Furthermore, within the trilin- 
ear ME interaction, even if two irreps are present, if they 
are in phase [i. e. if {4> — <^e ) /tt is an integer], then a 
spontaneous polarization does not arise^ When cooling 
from the paramagnetic phase into the (/, 0, /) phase, one 
expects only a single irrepi^ 3 - Upon further cooling, sys- 
tems that follow the scenario of Fig. la condense a sec- 
ond irrep and thereby 2 ^ induce ferroelectricity. When we 
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TABLE VIII: Character table for the point group for the 
125's. T a , where a = x,y,z are vector irreps. The next- 
to-last row gives the characters of the 34-dimensional re- 
ducible representation V and the last row gives those of the 
18-dimensional reducible representation, IV. 





E 


m bc 


m ac 


m ab 


1 


2a 


2„ 


2 C 




1 


1 


1 


1 


1 


1 


1 


1 


r* 


1 


-1 


1 


1 


-1 


1 


-1 


-1 


r, 


1 


1 


-1 


1 


-1 


-1 


1 


-1 


i+ 


1 


1 


1 


-1 


-1 


-1 


-1 


1 




1 


1 


-1 


-1 


1 


1 


-1 


-1 


r 

x xz 


1 


-1 


1 


-1 


1 


-1 


1 


-1 


r 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


r 

J- xyz 


1 


-1 


-1 


-1 


-1 


1 


1 


1 


r 


34 


4 


2 


4 


10 


4 


10 


4 




18 


6 


2 
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have both irreps of the wave vector present, their relative 
phase [0(r e ) — (f)(Y )}/-K is usually fixed by fourth order 
terms in the magnetic free energy to be nonintegral^^ in 
which case no choice of origin will simultaneously make 
both irreps inversion invariant. This situation is remi- 
niscent of TMOi or NVO^ and was previously noted in 
connection with second harmonic generation.— Finally, 
from Eq. (|53[) one sees that even when two irreps are 
present, if the order parameters of the two wave vectors 
q+ and q~ have the same magnitude, the spontaneous 
polarization could vanish. (This probably corresponds 
to the spirals of the two wave vectors having opposite 
helicity4£) 



B. Higher Order ME Coupling 

Sergienko et alM- have pointed out the existence of 
higher order terms in the ME coupling, in particular 
terms quartic in the order parameters. As they indicate, 
these terms have the potential to induce a spontaneous 
polarization in direction(s) different from those of the 
trilinear ME coupling. For the so-called 113 compounds 
(such as HoMnC>3, which they consider), these terms usu- 
ally do not come into play in view of the anisotropy of 
the terms in the purely magnetic free energy which are 
quartic in the order parameters. (See citation 28 of Ref. 
|24| .) Here the situation is different: the quartic order- 
parameter anisotropy is much more complicated for the 
125's, so that these higher order ME terms may come 
into play, although, as mentioned, their effect may be 
small. We start by first considering terms which strictly 
conserve wave vector. Later, we will investigate the corre- 
sponding Umklapp terms which only conserve wave vector 
to within a nonzero reciprocal lattice vector. 

To construct this ME interaction we need to construct 



quartic terms in the order parameters which transform 
like a vector. To avoid complications, it is simplest to 
use the following approach suggested by Mukamel^S The 
idea is to first find the number of such vector represen- 
tations by using the character tables to determine how 
many times each vector irrep is contained in the reducible 
representation formed by the basis functions of all fourth 
order terms. The 34 fourth order terms are the nine dis- 
tinct terms of the form 

cr fc (q + )CTi(q + )CT m (q + )V„(q + )* , (54) 
the nine distinct terms of the form 

cr fe (q_)o7(q_)o- m (q_)V n (q_)* , (55) 
and the 16 terms of the form 

CTfc(q+)cr;(q-)crm(q+)*o„(q_)* , (56) 

where k, I, m, and n assume the values o and e. The 
character table for the irreps of the point group of Pbam 
and that for the representation T generated by the quar- 
tic terms are given in Table IVIIII The characters of the 
representation T for each operator are obtained by tak- 
ing the trace of the operator in the 34 dimensional vector 
space under consideration. 

Then, we find the number of times n(T a ) that T a is 
contained in T is given by the scalar products of the char- 
acter vectors given in Table IVHll as 51 

n(r x ) = (34-4 + 2 + 4- 10 + 4- 10 -4)/8 = 2 , 
n(T y ) = (34 + 4-2 + 4- 10 -4 + 10- 4)/8 = 4 , 
n(T z ) = (34 + 4 + 2-4- 10-4- 10 + 4)/8 = 2 . 

(57) 

We find the two x-like functions to be 

i 22,22 22 22 

4>x,l = V 3 V 5 + V 4 U 6 - V 2 V S - UjWiJ 

= Mq + K(q + )*] 2 + Mq_K(q_)*] 2 

-[ ( x e (q_) ( 7 (q_)*] 2 - Mq + K(q + )*] 2 , 

(f>x,2 = V3V4V 5 V 6 - ViV 2 V 7 V 8 

= Oo(q+)cr (q_)(7 e (q + )V e (q_)* 

-<J e (q + )er e (q_)cr (q + )*cr (q_)* . (58) 

The above are easy to check, at least apart from the com- 
plex phase factors which always combine to give unity. To 
be invariant under m ac we must have an even number of 
"o" 's and an even number of "e" 's. Note that to be odd 
under X, the form must be odd under complex conjuga- 
tion. To be even under m a b the form must be even under 
interchange of q+ and q . 

We find the four y-like functions to be 

(f>y,l = VxV 3 vl + V 2 V4,vl - vlv 6 V$ - v\v 5 V 7 

= Mq+)| 2 Mq+)cr e (q+)* - cr (q + )V e (q + )] 

+ Mq_)| 2 Mq-)cr e (q-)* - a (q_)V e (q_)] 
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4>y,2 = vlv 5 V 7 + v\vqV% - V 2 Vivl - V X V 3 V 7 

= K(q+)| 2 Mq+)o- e (q+)* - cr (q+)V e (q+)] 
+ |c r e(q-)| 2 [cr (q-)cr e (q-)* - o- (q_)V e (q_)] 
<t>y,3 = ViV 5 [v4V 6 - v s v 2 ] + v 2 v 6 [v 3 v 5 - V1V7] 

= k e (q+)| 2 ko(q-)o' e (q_)* - <7„(q_)V e (q_)] 

+Mq-)l 2 Mq+)o- e (q+)* - cro(q+)*o- e (q + )] 

4>y,4, = V 4 V$[viV 7 - V 3 V 5 ] + V 3 V 7 [v 2 V$ - V4,V 6 ] 

= ko(q-)| 2 [c r e(q+)cr (q+)* - cr (q + )cr e (q + )*] 

+|o-o(q+)| 2 [o- e (q-)o- (q-)* - o- (q_)<7 e (q_)*] . 

(59) 

These can be checked similarly. To be odd under m ac the 
"e" ! s and the "o" 's must both appear an odd number of 
times. 

We find the two z-like functions to be 
1 22 22,22 22 

9z,l = «3 U 5 ~ v i v 6 + V 2 V 8 ~ v l v 7 

= Mq + )a e (q + )*] 2 - [cx (q_K(q_)*] 2 

+ [a e (q_K(q_)*] 2 - [a e (q + K(q + )*] 2 , 

4>z,2 = V 2 V 3 V 5 V S - ViV4V 6 V 7 

= CT e (q_)CT (q + )o- e (q+)*CT (q_)* 

-cr e (q_)*cr (q + )*cr e (q + )c7 (q_) . (60) 

These can be checked similarly. To be odd under TO a f, 
the form must be odd under interchange of q + and q_ . 
The ME interaction of order c 4 is written as 

^ME = ^ ] Cn,~f4'~f,nP'y , (61) 

where the c„ i7 are unknown coefficients. Now we discuss 

how V^ 4 g affects the ME phase diagrams. First of all, if 
there is only a single irrep, either an "e" or an "o" , then 
this interaction vanishes. So in the (1,0,1) phase, which 
has only a single irrep,— we still have no spontaneous 
polarization. As mentioned in the introduction to this 
section, this higher order ME interaction may be small 
and difficult to observe. 



C. Umklapp ME Interactions 

Now we consider Umklapp terms relevant to the phase 
in which q z = 1/4 but q x 7^ 1/2. Here the reducible 
representation Ty is generated by the nine terms of the 
form 

CT fe (q + )CTi(q + )CT m (q_)*CT„(q_)*(54 92 ,i , (62) 

and the nine terms of the form 

£rfc(q_)a-;(q_)o- m (q + )*CT n (q + )*<5 4 g a! ,i . (63) 

The characters for Tjj are given in Table |VllTJ Then, we 
find the number of times n(T a ) that r Q is contained in 



T v to be£i 

n(T x ) = (18-6 + 2 + -0 + -0 - 6)/8 = 1 , 
n(T y ) = (18 + 6-2 + 0-0-0 + 0-6)/8 = 2 , 
n(T x ) = (18 + 6 + 2 -0-0-0-0 + 6)/8 = 4 . 

(64) 

We find the x-like function to be 

/ 22,22 22 22 

ipx,i = v 3 v 6 + u 4 u 5 - v 2 v 7 - v x v s 

= Mq + H(q+)*] 2 + [a (q_) ( r e (q_)*] 2 

-[<7 e (q_)o- (q_)*] 2 - [a e (q+K(q + )*] 2 , 

(65) 

the two y-like functions to be 

Ipy, 1 = V!V 3 vl + V 2 V 4 vl - vjv 5 V 7 - vlv 6 V 8 

= o- e (q+H(q-)*Mq+)er e (q_)* 
-CT e (q+)cr (q_)*] + <j e (q_)o- e (q + )* 
x[o- (q_)CT e (q+)* - cr e (q_)cr (q + )*] 

1py,2 = vlv 6 V S + vjv 5 V 7 - V 2 ViV 2 - viv 3 vl 

= cr (q + )(7 (q_)*[c7- (q + )cre(q-)* 
-cr (q_)*cr e (q + )] + cr (q_)cr (q + )* 
x[cr (q_)o- e (q + )* - <j (q + )*o- e (q_)] , 

(66) 

and the four z-like functions to be 

, 2222 

1f>Z,l = VlV 6 - V 2 V 5 

= [a e (q + )<7 e (q_)*] 2 - [ ( x e (q_) ( 7 e (q + )*] 2 

^z.2 = V 3 V% - v\v 2 

= Mq+y o (q_)*] 2 - Mq_)a (q + )*] 2 
1pz,3 = VlV 3 V 6 V$ - v 2 v 4 v 5 v 7 

= o- e (q+)^o(q+)o- e (q-)*cr (q_)* 

-a e (q_ ) a D (q_ ) a e (q+ ) * a (q+ ) * 

< 22 22 22,22 

VzA = v 3 v 6 - v 2 v 7 - v±v % + v i v 5 

= [a (q + ) ( r e (q_)*] 2 - [a e (q_K(q + )*] 2 

+ [a e (q + K(q_)*] 2 -Mq-K(q + )*] 2 .(67) 

The transformation properties of the ip a ,n can be checked 
just as we did for the <j) a ,n- The Umklapp ME interaction 
of order cr 4 is written as 

4e,u = ^lEC^."^' ( 68 ) 

where the cL „ are unknown coefficients. 

Clearly this interaction is only operative when q z is 
locked to the CM value q z = 1/4. This is therefore a gen- 
eralization of the term introduced by Betouras et al.^ 
but here we give the first analysis of the symmetry of this 
interaction. It is interesting to note that this interaction 
can induce a spontaneous polarization along the z-axis 
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even when only a single irrep is present. (Inspection of 
ip x i and indicates that this requires simultaneous 
condensation of order at wave vectors q±.) However, as 
mentioned in the introduction to this section, this higher 
order ME interactions may be small and difficult to ob- 
serve. 



COMPATIBILITY RELATIONS 



PARAM AGNET 



P-Pr 




FIG. 2: (Color online) Phase diagram (simplified from Ref. 
|2~I ) as a function of V and T for fixed q z near the multicriti- 
cal point. The phase with q x = 1/2 exists within a parabolic 
"tongue" whose apex is the multicritical point M, where the 
F e and F irreps interchange stability. The compatibility rela- 
tions we obtain apply in the vicinity of the multicritical point 
M. 



A first step to constructing a generic phase diagram 
for the 125's2^ is to understand how the wave functions 
behave near the phase transition between the phase with 
q x 1/2 and that for which q x = 1/2. In Fig. [2] we 
show a simplified version of this phase diagram for fixed 
q z . (However, to compare with experiment, the diagram 
for fixed q x is more relevant To avoid confusion we 
introduce a control parameter V such that when V = V C 
the wave vector which minimizes the inverse susceptibil- 
ity near the ordering transition has q x = 1/2, but when V 
deviates slightly from this critical value the ^-component 
of the selected wave vector is not exactly equal to 1/2. 
(We refer to the point V = V c as "the multicritical point" 
because to reach this point requires not only fixing the 
temperature to be at the ordering transition, but also, 
as shown in Fig. [21 one must fix V — V c by varying 
some other parameter, such as the pressure.) As we have 
seen, as the temperature is lowered into the ordered phase 
when V ^ V c , one of the ID irreps r e or r o at q x ^ 1/2 
condenses, whereas exactly at V — V c one condenses 
into a phase with q x = 1/2 which has only a single 2D 
irrepj*i2£ Accordingly, we now study the compatibility re- 
lation which must relate the wave functions of these two 
phases in the limit as we approach the multicritical point 
M for which V = V c . Experimentally, the phase transi- 
tion between the phase with q x = 1/2 and that having 



q x 1/2 has only been observed for q z = 1/4. How- 
ever, since the symmetry of the phases for q z — 1/4 is 
not different from that for q z ^ 1/4, we will leave q z as 
a free parameter which we consider to be incommensu- 
rate. Although the actual phase transition between the 
q z = 1/2 phase and the phase with q z 7^ 1/2 must be dis- 
continuous, the discontinuity vanishes in the limit when 
the multicritical point M in Fig. [5] is approached. In this 
limit, one may consider the transition to be continuous, 
and therefore it must be possible to express each basis 
function of the two irreps of the q x ^ 1/2 phase as a 
linear combination of the basis functions of the 2D irrep 
of the phase having q x = 1/2. We do this explicitly in 
order to find the relation between the order parameters 
of the two phases. This relation will be perturbatively 
modified as one goes deeper into the ordered phase. 



A. Wavefunctions near the Multicritical Point 

In Table IIXI we record the wave functions allowed by 
symmetry for the q x = 1/2 state, based on Table XVI 
of Ref. 9, which are modified in several ways. First of 
all, one has to include the corrections to the wave func- 
tions on sublattices 9-12, as described in an erratum^ 
Secondly, we translate all sites by (0, 0, 1/2). (This oper- 
ation has no effect because the induced change of phase 
can be absorbed into the order parameters.) Thirdly, we 
renumber the sublattices to make their positions equal to 
their counterparts in TableHIlto within a lattice constant. 
The final step was to translate sublattices through an in- 
teger number of lattice constants, as necessary, in order 
to bring them back into the unit cell. In this last opera- 
tion sublattice n was translated through A„, where Ai = 
(0,1,0), A ± = (1,0,0), A 9 = A10 = An = (0,0,1), and 
A 12 = (1,1,1). The result of this operation was to in- 
troduce a multiplicative factor X n — exp[27riq • A ra ] to 
all components of the nth sublattice. Thereby we obtain 
the results shown in Table 11X1 

Near the multicritical point M the critical spin wave 
function v I / 9x= i/2 (for a fixed value of q z and q x — 1/2) is 
a linear combination of a\ times the basis functions of the 
first column of Table IIXI plus a 2 times the basis function 
of the second column of Table HXl Alternatively, near the 
multicritical point M for q x ^ 1/2 phase, this spin wave 
function can be formed within the space in which the two 
ID irreps, T e and r o , are considered degenerate for the 
fixed value of q z . In this limit the wave function ^ qx7 ti/2 
of the ID irrep phase is given by a linear combination 
of the basis functions associated with the four order pa- 
rameters = lim^o c s [±(l/2 — S), 0, q z ], where s is e 
or o. These basis functions are given in Tables IIVI and 
IVl Equating ^^=1/2 and ^ fe ^i/2 gives, with, as before, 

= (-1, 1, -1) and p a = (-1, -1, 1), 

01 r la + o- 2 r 2a 

= c ( ts„,i + (7+u Qj i + er~ p a s* al + a~ p a U al (69) 
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TABLE IX: Normalized spin functions (i. e. Fourier coef- 
ficients) within the unit cell of e. g. TbM^Os for wave 
vector (i,0, g). Here the r nx , r ny and ir nz are real, the z's 
are complex, and A = exp(27rig z ) where q z is in rlu's. The 
x, y, and z components of each Fourier vector are listed in 
the corresponding box. The actual spin structure is a linear 
combination, o\ times the first column plus cri times the sec- 
ond column, where the cr's are complex order parameters and 
the entries in each column are normalized so that the sum of 
their absolute squares is unity. 



Spin 


cri 


CTl 


Spin 


cri 


cri 




Tlx 


Tlx 




Zx 


Zx 


S(q,l) 


n y 


r 2y 


S(q,7) 


Zy 


Zy 




ru 


r 2 z 




-Z z 


Zz 




rix 


—ria; 




Zx 


— Z x 


S(q,2) 


-r 2y 


ri y 


S(q,8) 


— Zy 


Zy 




riz 


—riz 




Zz 


Zz 




rix 


—rix 




r.xA 1 ' 2 


r.xA 1 ' 2 


S(q,3) 


-r ly 


ri y 


S(q,9) 


r 5y A- 1/2 


r 6y A^ 2 




—r\ z 


riz 




r 5 zA 1/2 


rezA 1 ? 2 




-r 2x 


-rix 




rexA 1 ? 2 


—r 5x A 1/2 


S(q,4) 


-r 2y 


-riy 


S(q, 10) 


-r 6y A^ 2 


r 5yA 1/2 




riz 


riz 




rezA^ 2 


-r.zA 1 ' 2 




Zx 


Zx 




rsxA^' 2 


-r 6x A 1/2 


S(q,5) 




z* 

6y 


S(q,ll) 


—r 5y A 1/2 


r 6y A" 2 




-z% 


-Z*z 




-r 5z A^ 2 


rezA 1 ' 2 




Zx 


z x 




-r,xA^ 2 


-rsxA 1 ' 2 


S(q,6) 


Zy 


Zy 


S(q, 12) 


-r, y /A" 2 


-r 5y A^ 2 




Z*z 






rez/A 1 ' 2 


rszA 1 ' 2 



-Ca[crir 2 Q - tT 2 r la ] 

= {-i£,a)vi ~s a ,i + (iia)cri ~u a ,i 

) a e s a 1 + ( _ ^aPa)c o u *a 1 (70) 

= (*£a)o-+s*,i + (-if a )<7+ u* aA 

+ {~iPa£,a)<y~ S Q ,1 + {i£, a p a )<T~ U a ^i (71) 

= ct S a,l + a o u * a ,l + °e PaS a ,l + p a U a .i (72) 
£, a Z* a (p a (7l + CT 2 ) 

= vtSa,2 + &o u a,2 + &e Pa s a,2 + &o Pa u a,2 (73) 
z a( a l - Pa 0"2 ) 

Pa Sql,2 + (^a)o'o PaUa,2 (74) 

z a (-p a cri + (T 2 ) 
= iZ, a ot s * a ,2 + {-^a)(JoK.2 
+ {-i£ a )(T~ p a S* a 2 + (i£, a )<7o PaK,2 (75) 



+ PaO- 2 ) 

= <rt s *a,2 + + °ePaS* a . 2 + <?o PaK,2 (76) 

(o-ir- 5iQ! + a 2 r e ^)A 1/2 

= ^s a ,3 + cr+n Q , 3 + <r~Ap a s* a 3 + a~ Ap Q u* )£ (77) 
-£a(0ir 6 , a - (T 2 r 5: a )A 1/2 
-\-{^a)Ap a (T~ s* a3 + {-i^ a )Kp a a~u* a<i (78) 

£aPa(crir 5 , a - 0'2?'6,q)A 1/2 

= A(t6>+«* )3 + A(-i^ a )a+u* a3 

+ (-i£a)Pa<re s a,3 + (*fa) Pa<J~ U Q ,3 (79) 

f J (\ 

)A l/2 

= Acr+s* 3 + Acr+u* 3 + cr~p a s at3 + <J~p a u at3 (fiO) 

where erf = lim^o c s [±(l/2 — S),Q,q z ], where s is e or 
o. 



% \Qx) 




Qx 



1/2 

FIG. 3: (Color online) The two lowest eigenvalues of X~ 1 {Qx) 
which are degenerate for V = V c - 

B. Symmetry of the Multicritical Point 

From the above equations we expect to obtain a re- 
lation between the order parameters of the phase with 
q x = 1/2 and that with q x ^ 1/2 arbitrarily close to the 
multicritical point M. Presumably, giving the values of 
<7g and <7q will determine the values of cri and 02, but 
having the values of <J\ and a 2 we can not expect to de- 
termine the four parameters erf- and . Accordingly, 
we now study the basis functions for T e and r o and show 
that they are related in the limit when q x — > 1/2. To see 
this we will analyze the behavior of the inverse suscepti- 
bility as a function of Q Xl the x-component of the wave 
vector when the temperature is just above the tempera- 
ture at which magnetic order appears and for V close to 
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the critical value V c at which the minimum of the inverse 
susceptibility as a function of Q x occurs for Q x = 1/2. 
Note that the inverse susceptibility has 36 branches, each 
one corresponding to an eigenvalue of the inverse suscep- 
tibility matrix. Here we need consider only the two lowest 
branches of the inverse susceptibility. These lowest two 
eigenvalues arise out of a two by two submatrix which we 
now analyze for Q x = 1/2 + k x and V = V c + y for small 
k x and y. For y = this submatrix is of the form 



a(T- 







bk 2 x 

a(T - T c ) + bk 2 x 



(81) 



where k x = Q x — 1/2, a and b are constants, and T c is 
the temperature at which order first develops. Here and 
below we work only to order k x . This form is dictated by 
the fact that the inverse susceptibility has to be two-fold 
degenerate, have its minima at k x = 0, and the spectrum 
has to be independent of the sign of k x (in view of the 
existence of the symmetry element m/, c ). Thus the two 
lowest branches in the eigenvalue spectrum of the inverse 
susceptibility as a function of Q x are as shown in Fig. [3] 



O (1/2-8 ,0,q z ) 
(u f u 2 ,u 3 ) 



O e (1/2-5 ,0,q 2 
(s , s , s ) 

v 1' 2' 3 ' 




O (-1/2+5 ,0, qz ) 
(pU*,pU 2 ,pAl|) 



a e (-1/2+5 ,0,q z ) 
(PS*, p S 2 ,pAS*) 



1/2 



Qx 



FIG. 4: (Color online) The two lowest branches of eigenvalues 
of X _1 (9a0 with their wave functions indicated. Note that the 
labels e and o refer to the eigenvalues rather than the branch 
of the spectrum. We assume that the wave functions at q + for 
irreps F e and r o are given in terms of s n and u„ , respectively, 
as listed in Tables ITVl and IVl respectively. Here s n = s a , n , 
u n = u a ,n> and p = p a . Then the wave functions at q — q~ 
are obtained in terms of those at q = q + according to Table 
IVII In general, the o and e wave functions are unrelated. 
However, as V — > V c , q + — q~ — ► and the two parabolas come 
into coincidence. In this situation the points corresponding 
to <r e and o come into coincidence. Therefore by continuity 
on either the right-hand or the left-hand parabola the a and 
<r e wave function on the same parabola become equal, leading 
to Eq. 



Next consider allowed terms which are linear in y but 
have an unspecified dependence on k x . These will give 



X 



bk 2 x +c(k x )y 
d(k x )*y t 



d(k x )y 
bkl + e{k x )y 



,(82) 



where c{k x ) and e(k x ) are real and r = a(T—T c ). For the 
spectrum to be the same for both signs of k x , c(k x )+e(k x ) 



must be an even function of k x . The term in [c(k x ) + 
e{k x )\ independent of k x leads to an allowed dependence 
of T c on y and the term of order k x leads to an allowed 
dependence of the coefficient b on y, so, in effect, up to 
order kl we have 



- bk 2 x + c'yk x 
d(k x )*y 



d(k x )y 
bkl - c'yk x 



(83) 



where now r and b have an allowed, but unimportant, 
dependence on y. Now consider the dependence of d{k x ) 
on k x . Suppose that d{k x ) were nonzero for k x = 0. This 
would imply that the minimum in the inverse suscepti- 
bility occurred for k x = 0, but that the eigenvalues were 
not degenerate. This contradicts group theory. So the 
generic case is that d(k x ) — (3k x + 0(k x ). Then the two 
eigenvalues are 



X± = t + bk x ± yk 3 



(84) 



This leads to two parabolic branches of the inverse sus- 
ceptibility with minima symmetrically displaced away 
from Q x = 1/2 by an amount linear in V — V c , as shown 
in Fig. [U As shown there, the left parabola at Q x = q + 
is associated with r e and is parametrized by the s's and 
the right parabola at Q x = q + is associated with r o and 
is parametrized by the it's. The corresponding basis func- 
tions are given explicitly in Tables IIVI and [Vj But the 
basis functions for r o and T e at q~ are related, respec- 
tively, to r o and T e at q + according to Table PVTl and this 
is indicated in Fig. [3J These eigenfunctions of the inverse 
susceptibility depend on wave vector, of course. But as 
V — > Pc, the two parabolas come into coincidence with 
their minimum at Q x = 1/2, and the points governed by 
cr e and a Q on the same parabola approach one another. 
Then in this limit, by continuity on the same parabola 
we obtain 

U a ,l = PaS* l, U Qi2 = PqS Qj2 , Ua,3 = PqAs* 3 . (85) 

It should be remarked, that this multicritical point is 
not a Lifshitz point, — At a Lifshitz point the coefficient 
of k^ in the inverse susceptibility vanishes. Here, in the 
generic case, this coefficient is nonzero, but the coefficient 
of k x , which here is allowed because of the double degen- 
eracy, vanishes. Furthermore, the Lifshitz point sepa- 
rates a regime of CM order from that of IC order. Here 
CM order (at the paramagnetic phase boundary) only 
occurs at a point (where the coefficient of k x changes 
sign thereby exchanging the instabilities of the two ID 
irreps). 



C. Compatibility Equations 

Using the relation between the u's and the s's, we see 
that Eqs. dMD-dSOl) become 



ctiTiq + a 2 r 2a = & s ai i + a p a s a l 



(86) 
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ci^2q — & 2 r\u = icr + s Qi i — ip a a s* j , (87) The order parameters are related by 



0\T\ a - a 2 T 2a = ip a <T + S* a l - %0 S a ,l 



(89) 



2/0(0-1 + PaCT 2 ) = (/5 Q cr + +cr )s a ,2 , (90) 



Vai^l ~ Pa<J 2 ) = - ip„(7 )s Q ,,2 , (91) 



- p a <7 2 ) = (-«p Q cr + +icr )s* 2 , (92) 



y a {<ri + Pacr 2 ) = (-er+-cr p a )s* c 



(93) 



uir^. a + cr 2 r 6iQ = er + s Qi3 A 1/2 

+ ( t-A 1 / 2 PqS ^ 3 , (94) 

cnre,a - o- 2 r 5iQ = ia + s Qj3 A~ 1/2 

-iAVa paff - a * j3 _ (g5) 

oir 5 , a - (T 2 r 6 , a = ip a cr + s* a j3 A 1/2 

-z^-^.gA- 1 / 2 , (96) 



crir 6 ,a + o- 2 r 5 . a = p a cr + s* 3 A 1/2 

+ C r-,S Qi3 A- 1 / 2 



(97) 



where y a — £, a z a and = a^+aT . These equations are 
strongly overdetermined. Accordingly, the fact that they 
have a solution is evidence that the wave functions which 
formed the input to this calculation are correct. (Indeed, 
in order to arrive at a solution, it was necessary to correct 
an error in the table of wave functions of Ref. 9.) These 
equations have the solution for the wave functions of the 
2D irrep phase in terms of those of the ID irrep phase as 



n, Q = [e™ /4 s a ,i - Pa e- l ^ 4 s* a J/V2 



(98) 



r 2 , a 



-e-™/ 4 s aA +p a e™/ 4 s* al ]/V2 , (99) 



-e OT/4 + p a e- M / 4 ]4 



(100) 



r 5 , a = [e^A-^Vs " Pae-^A^V 3 }/V2 ,(101) 



r 6 , a = [-e-^A-VVg + Pae^A 1 /^* 3 ]/y/2 (102) 



e-^ 4 a 2 ]/V2 



<j- = [-e-^/Vx 



The inverse transformation is 



e™/ 4 a 2 ]/V2 



(103) 



= \e-™ /A a+-e l 



g 2 = [-e l7r ^a+ 



■/V-]/V2 

M / 4 a-]/\/2 



(104) 

A strong check on these results is that the r nx and r ny 
are real (p x = p y = — 1) and r nz is imaginary (p 2 = 1), all 
as required by the symmetry analysis of the CM phased 

These results show how the order parameters of the 
2D irrep are related to the order parameters of the ID 
irreps. One should also note that by continuity, if the IC 
phase has a spontaneous polarization as q z — > 1/4, the 
CM phase should also have one, and vice versa. This is 
ensured by the fact that 

|<ti| 2 -M 2 = i[(at + a-)(a-+a+T 

a- )>-+<!+)] . (105) 

Now we only keep terms which conserve wave vector 
when we go away from q x = 1/2, in which case 



\°2\ 



-cr„ cr„ 



(106) 



Thus the ME interaction of Eq. ([53]) goes smoothly 
into the ME interaction in the CM state&^Ml 



= r[M 2 - H 2 ]P 6 



VI. CONCLUSION 



(107) 



We have performed a representation analysis of the 
magnetic order for the IC phase of the RMn 2 Os series by 
including inversion symmetry, thereby reducing by about 
half the number of degrees of freedom allowed for mag- 
netic ordering. Our results emphasize that a full inclu- 
sion of inversion symmetry is necessary to determine the 
magnetic structure and associated order parameters, not 
only in multiferroics, but also in a wide range of mag- 
netic materials. We have also determined the physically 
important order parameters and have analyzed the trans- 
formation properties which they inherit from the wave 
functions. Using these symmetry properties we have ana- 
lyzed the magnetoelectric interaction responsible for the 
simultaneous magnetic and dielectric phase transitions. 
The lowest order magnetoelectric interaction, which is 
bilinear in the magnetic order parameters, explains the 
observed direction of the spontaneous polarization. We 
have shown that higher order and Umklapp magnetoelec- 
tric interactions (which are quartic in the spin variables) 
can induce nonzero values for all components of the spon- 
taneous polarization, but since the order parameters are 
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small in the relevant phases and since microscopic mech- 
anisms tend to involve terms quadratic in the spin vari- 
ables, these anomalous components to the spontaneous 
polarization may be very difficult to observe. We have 
also explicitly obtained the compatibility relations for the 
transition between the IC phase and the CM phase (or 



more generally the phase where the x-component of wave 
vector is locked to its CM value). 
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